Entanglement entropy at infinite randomness fixed points in higher dimensions 
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The entanglement entropy of the two-dimensional random transverse Ising model is studied with 
a numerical implementation of the strong disorder renormalization group. The asymptotic behavior 
of the entropy per surface area diverges at, and only at, the quantum phase transition that is 
governed by an infinite randomness fixed point. Here we identify a double-logarithmic multiplicative 
correction to the area law for the entanglement entropy. This contrasts with the pure area law valid 
at the infinite randomness fixed point in the diluted transverse Ising model in higher dimensions. 

PACS numbers: Valid PACS appear here 



Extensive studies have been devoted recently to under- 
stand ground state entanglement in quantum many-body 
systems []]]. In particular, the behavior of various entan- 
glement measures at /near quantum phase transitions has 
been of special interest. One of the widely used entan- 
glement measures is the von Neumann entropy, which 
quantifies entanglement of a pure quantum state in a bi- 
partite system. Critical ground states in one dimension 
(ID) are known to have entanglement entropy that di- 
verges logarithmically in the subsystem size with a uni- 
versal coefficient determined by the central charge of the 
associated conformal field theory Q • Away from the crit- 
ical point, the entanglement entropy saturates to a finite 
value, which is related to the finite correlation length. 

In higher dimensions, the scaling behavior of the entan- 
glement entropy is far less clear. A standard expectation 
is that non-critical entanglement entropy scales as the 
area of the boundary between the subsystems, known as 
the " area law" 0, [j] . This area-relationship is known to 
be violated for gapless fermionic systems [5j in which a 
logarithmic multiplicative correction is found. One might 
suspect that whether the area law holds or not depends 
on whether the correlation length is finite or diverges. 
However, it has turned out that the situation is more 
complex: numerical findings and a recent analytical 
study || have shown that the area law holds even for 
critical bosonic systems, despite a divergent correlation 
length. This indicates that the length scale associated 
with entanglement may differ from the correlation length. 
Another ongoing research activity for entanglement in 
higher spatial dimensions is to understand topological 
contributions to the entanglement entropy 

The nature of quantum phase transitions with 
quenched randomness is in many systems quite different 
from the pure case. For instance, in a class of systems 
the critical behavior is governed by a so-called infinite- 
randomness fixed point (IRFP), at which the energy scale 
e and the length scale L are related as: lne ~ with 
< V < 1- I n these systems the off-critical regions 
are also gapless and the excitation energies in these so- 



called Griffiths phases scale as e ~ L~ z with a nonuni- 
versal dynamical exponent z < oo. Even so, certain ran- 
dom critical points in ID are shown to have logarith- 
mic divergences of entanglement entropy with universal 
coefficients, as in the pure case; these include infinite- 
randomness fixed points in the random-singlet universal- 
ity clas s fl2 . 13L 14, 15, 1|| and a class of aperiodic singlet 
phases (171 ]. 

In this paper we consider the random quantum Ising 
model in two dimensions (2D), and examine the disorder- 
averaged entanglement entropy. The critical behavior of 
this system is governed by an IRFP 0, llll | implying that 
the disorder strength grows without limit as the system is 
coarse grained in the renormalization group (RG) sense. 
In our study, the ground state of the system and the 
entanglement entropy are numerically calculated using a 
strong-disorder RG method [HI, [ll| , which yields asymp- 
totically exact results at an IRFP. To our knowledge this 
is the first study of entanglement in higher dimensional 
interacting quantum systems with disorder. 

The random transverse Ising model is defined by the 
Hamiltonian 



H 



E J iJ a i a j 

(id) 



(1) 



Here the {erf } are spin- 1/2 Pauli matrices at site i of an 
L x L square lattice with periodic boundary conditions. 
The nearest neighbor bonds Jy(> 0) are independent 
random variables, while the transverse fields hi(> 0) are 
random or constant. For a given realization of random- 
ness we consider a square block A of linear size £, and 
calculate the entanglement between A and the rest of 
the system B, which is quantified by the von Neumann 
entropy of the reduced density matrix for either subsys- 
tems: 

S = -Tr(p A log 2 p A ) = -Tr(p B log 2 p B ). (2) 

The basic idea of the strong disorder RG (SDRG) ap- 
proach is as follows IS, HI: The ground state of the sys- 
tem is calculated by successively eliminating the largest 
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FIG. 1: (color online). An example of typical ground state 
in the random quantum Ising model (a) in ID, and (b) in 2D; 
it contains a collection of spin clusters of various sizes, which 
are formed and decimated during the RG. The entanglement 
of a block (shaded area) is give by the number of decimated 
clusters (indicated by red loops) that connect the block with 
the rest of the system. 



local terms in the Hamiltonian and by generating a new 
effective Hamiltonian in the frame of the perturbation 
theory. If the strongest bond is , the two spins at i and 
j are combined into a ferromagnetic cluster with an effec- 
tive transverse field h^j) — -j- 2 -- If, on the other hand, 
the largest term is the field hi, the spin at i is decimated 
and an effective bond is generated between its neighbor- 
ing sites, say j and k, with strength Jjk — J 'l .** . After 
decimating all degrees of freedom, we obtain the ground 



state of the system, consisting of a collection of indepen- 
dent ferromagnetic clusters of various sizes; each cluster 
of n spins is frozen in an entangled state of the form: 



(3) 



In this representation, the entanglement entropy of a 
block is given by the number of clusters that connect 
sites inside to sites outside the block [Fig. [T] . We note 
that correlations between remote sites also contribute to 
the entropy due to long-range effective bonds generated 
under renormalization. 

In ID the RG calculation can be carried out analyti- 
cally and the disorder-averaged entropy S# of a segment 
of length £ has been obtained as Si = log 2 I [13] • 
In higher dimensions the RG method can only be imple- 
mented numerically. The major complication in this case 
is that the model is not self-dual and thus the location of 
the critical point is not exactly known. To locate the crit- 
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FIG. 2: (color online). The distribution of the last deci- 
mated effective log-fields In h^, , and the distribution of the 
last decimated effective log-bonds In Joo in the RG calcula- 
tions. At ho = 1.175, the distributions, shown in (a) and (b), 
get broader with increasing system sizes, indicating the RG 
flow towards infinite randomness, i.e. the system is critical. 
A scaling plot of the data in (a) using energy-length scaling 
Inhoc ~ with tp = 0.55 is presented in (c). The solid line 
is just a guide to the eye. The subfigures (d) and (e) show 
the log-field distribution at ho = 1.18 and the log-bond dis- 
tribution at ho = 1.17, respectively; the distributions show 
a power-law decaying tail in the low energy region, which is 
clear evidence that the system is in the Griffiths phases. 



ical point, we can make use of the fact that the excitation 
energy of the system has the scaling behavior In e ~ 
at criticality, while it follows e ~ L~ z in the off-critical 
regions. In the numerical implementation of the SDRG 
method, the low-energy excitations of a given sample can 
be identified with the effective transverse field /i M of the 
last decimated spin cluster, or with the effective coupling 
Joo of the last decimated cluster-pair. 

In our implementation we set for convenience the 
transverse fields to be a constant ho and the random 
bond variables were taken from a rectangular distribu- 
tion centered at J = 1 with a width A = 0.5. The critical 
point was approached by varying the single control pa- 
rameter /iq. Although this initial disorder appears to be 
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FIG. 3: (color online). Left panel: The disorder averaged 
block entropy per surface unit Se/£ vs. the linear size of the 
block £ for a system size L — 64 for various values of ho- We 
observe that the entropy for £ = L/2 reaches its maximum 
at the critical point h c — 1.175 (cf. Fig [2]). Right panel: 
The block entropy per surface area vs. ln£ on a log-scale for 
different system sizes L at the critical point. The data show 
a straight line (guided by the dashed line), corresponding to 
the scaling obeying the area law with a double-logarithmic 
correction, as given in Eq. Q. 



weak, the renormalized field and bond distributions be- 
come extremely broad even on a logarithmic scale [Fig. [5] 
at the critical point ho = h c = 1.175. This indicates the 
RG flow towards infinite randomness. Slightly away from 
the critical point, both in the disordered Griffiths-phase 
with ho = 1.18 and in the ordered Griffiths- phase with 
ho = 1.17, the distributions have a finite width and obey 
quantum-Griffiths scaling h x ~ L~ z . At the critical 
point one has IRFP scaling In hoo ~ Lr and we estimate 
the scaling exponent as ip = 0.55, quite close to the value 
tp = 0.5 for the ID case 

Now we consider the entanglement entropy near the in- 
finite randomness critical point. To obtain the disorder- 
averaged entanglement entropy Sg of a square block of 
size £, we averaged the entropies over blocks in different 
positions of the whole system for a given disorder real- 
ization and then averaged over a few thousand samples. 
In Fig. [3] we show the entropy per surface unit Sg/£ = st 
for different values of ho. This average entropy density 
is found to be saturated outside the critical point, which 
corresponds to the area law. At the critical point sg in- 
creases monotonously with £, and the numerical data are 
consistent with a log-log dependence: 



S e ~ £log 2 log 2 £ 



(4) 



as illustrated in Fig. [3] In this way we have identified an 
alternative route to locate the infinite randomness criti- 
cal point: it is given by the field ho for which the average 
block entropy at I = L/2 is maximal. Indeed the nu- 
merical results in Fig. [3] predict the same value of h c as 
obtained from the scaling of the gaps. We note that the 



same quantity, the position of the maxima of the average 
entropy, can be used for the random quan tum Ising chain 
to locate finite-size transition points [21] . 

The log-log size dependence of the average entropy in 
Eq. (HJ) at criticality is completely new; it differs from the 
scaling behavior observed in 2D pure systems, like the 
area law, Sg ~ £, for critical bosonic systems HI, or 
a logarithmic multiplicative correction to the area law, 
St ~ I log 2 £, as found in free fermions 

BUSH- This 

double-logarithmic correction can be understood via a 
SDRG argument: In the ID case a characteristic length 
scale r at a given RG step is identified with the aver- 
age length of the effective bonds, i.e. the average size 
of the effective clusters. At the scale r(< £) the frac- 
tion of the total number of spins, rv, that have not been 
decimated is given by n r ~ 1/r 18]; these active (i.e., 
undecimated) spins have a finite probability to form a 
cluster across the boundary of the block (a segment £ in 
the ID case) and thus to give contributions to the en- 
tanglement entropy. Repeating the renormalization until 
the scale r ~ £, the contributions to the entropy are 
summed up: St ~ / drn r ~ ln£, leading to the log- 
arithmic dependence of the ID model [l2j. For the 2D 
case with the same type of RG transformation with a 
length scale r < £, the fraction of active spins in the 
renormalized surface layer of the block is n r ~ l/r. Here 
we have to consider the situation in which some of these 
active surface spins would form clusters within the sur- 
face layer and thus contribute zero entanglement entropy; 
the number of the active spins that are already engaged 
in clusters on the surface at RG scale r is proportional 
to lnr, as known from the ID case, and only 0(1) of 
the active surface spins would form clusters connecting 
the block with the rest of the system. Consequently, 
the entropy contribution in 2D can be estimated as: 
St ~ J rg dV n r / In r ~ £ In In £, i.e. a double-logarithmic £- 
dependence, as reflected by the numerical data in Fig. [3] 

Based on the SDRG argument described above, the 
double-logarithmic correction to the area law appears to 
be applicable for a broad class of critical points in 2D with 
infinite randomness. For instance, the critical points of 
quantum Ising spin glasses are believed to belong to the 
same universality class as ferromagnets since the frustra- 
tion becomes irrelevant under RG transformation, and 
the same type of cluster formations as observed in our 
numerics for the ferromagnet is expected to be generated 
during the action of the RG. The entanglement entropy 
at the IRFP is completely determined by the cluster ge- 
ometries occurring during the SDRG. 

Another type of IRFP in higher dimensions occurs 
in the bond-diluted quantum Ising ferromagnet: The 
Hamiltonian is again given by (P), but now Jij = with 
probability p and Jij = J > with probability 1 — p. At 
percolation threshold p = p c there is a quantum critical 
line along small nonzero transverse fields, which is con- 
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FIG. 4: (color online). The entropy per surface area 
Si/ 1 = S£ vs. £ near the percolation threshold p c — 0.5 for 
the 2D bond-diluted Ising model at small transverse fields for 
L — 512. The curves converge to finite values for £ — > oo, 
corresponding to the area law. The inset shows ~§t — ~s oc as 
a function of £. Soo is estimated from s L ^ 2 at L = 512. The 
dashed line corresponds to £~ . 



trolled by the classical percolation fixed point, and the 
energy scaling across this transition line obeys hie ~ L^, 
implying an IRFP |20(. The ground state of the system 
is given by a set of ordered clusters in the same geom- 
etry as in the classical percolation model - only nearest 
neighboring sites are combined into a cluster. In this 
cluster structure, the block entropy, determined by the 
number of the clusters connecting the block and the rest 
of the system, is bounded by the area of the block, i.e. 
Si ~ £ d ^ 1 with d being the dimensionality of the sys- 
tem. To examine this, we determined the entanglement 
entropy by analyzing the cluster geometry of the bond- 
diluted transverse Ising model. Fig. 0] shows our results 
for the square lattice, which follow a pure area-law with 
an additive constant: St = al + b + 0(1/ 1). 

To summarize, we have found that the entanglement 
properties at quantum phase transitions of disordered 
systems in dimensions larger than one can behave quite 
differently. Generalizing our arguments for the 2D case, 
we expect for the random bond transverse Ising systems 
a multiplicative d-fold logarithmic correction to the area 
law in d dimensions at the critical point, whereas for 
diluted Ising model at small transverse fields the area 
law will hold in any dimension d > 1 at the percolation 
threshold. Although both critical points are described 
by infinite randomness fixed points, the structure of the 
strongly coupled clusters in both cases is fundamentally 
different, reflecting the different degrees of quantum me- 
chanical entanglement in the ground state of the two sys- 
tems. This behavior appears to be in contrast to one- 
dimensional systems governed by IRFPs [12J . 

Other disordered quantum systems in higher dimen- 
sions might also display interesting entanglement prop- 



erties: For instance, the numerical SDRG has also been 
applied to higher dimensional random Heisenberg anti- 
ferromagnets which do not display an IRFP 22]. The 
ground states involve both singlet spins and clusters with 
larger moments; therefore, we expect the correction to 
the area law to be weaker than a multiplicative loga- 
rithm and different from the valence bond entanglement 
entropy in the Neel Phase [23j |. 
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